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ABSTRACT 

Forbes et al. recently used the Hubble Space Telescope to localize hundreds of 
candidate star clusters in NGC1023, an early-type galaxy at a distance of 11.1 
Mpc. Old stars dominate the light of 92% of the clusters and intermediate-age 
stars dominate the light of the remaining 8%. Theory predicts that clusters with 
such ages can host intermediate-mass black holes (IMBHs) with masses Mbh 
10 5 M 0 . To investigate this prediction, we used 264 s of 5.5 GHz data from the 
Karl G. Jansky Very Large Array (VLA) to search for the radiative signatures of 
IMBH accretion from 337 candidate clusters in an image spanning 492" (26 kpc) 
with a resolution of O'MO (22 pc). None of the individual clusters are detected, 
nor are weighted-mean image stacks of the 311 old clusters, the 26 intermediate- 
age clusters, and the 20 clusters with stellar masses M* > 7.5 x 10 5 M 0 . The 
clusters thus lack radio analogs of HLX-1, a strong IMBH candidate in a cluster 
in the early-type galaxy ESO 243-49. This suggests that HLX-1 is accreting gas 
related to its cluster’s light-dominating young stars. Alternatively, the HLX-1 
phenomenon could be so rare that no radio analog is expected in NGC1023. 
Also, using a formalism heretofore applied to star clusters in the Milky Way, 
the radio-luminosity upper limit for the massive-cluster stack corresponds to a 
3 a IMBH mass of Mbh (massive) < 2.3 x 10 5 M 0 , suggesting black-hole mass 
fractions of MBH(massive)/M* < 0-05 — 0.29. 
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Subject headings: black hole physics — galaxies: individual (NGC 1023) — galax¬ 
ies: star clusters: individual (NGC 1023, ESO 243-49 HLX-1) — radio continuum: 
general 


Motivation 


Constraints on populations of low-mass black holes (BHs) are important for understand¬ 
ing formation channels for seed BHs in the early Uni verse a nd the demographics of relic seed 
BHs in the local Universe (reviewed by Greene 2012 : Volonterill2012[ l. Such seeds are referred 
to as intermediate-mass black holes (IMBHs) because they occupy the gap in mass between 
the well-studied stellar-mass BHs with M B h ~ 10 M e and the well-established supermassive 
BHs with Mbh > 10 6 M©. In addition, if the scaling relations between central BHs and their 
host stellar systems differ between the intermediate-mass and supermassive ends, then this 


could provide clues about the co-evoluti on, or not, of these entities (e.g., Kiao et al.l 1201 1 
Jiang et al. 2011 : Kormendv fe Ho 2013 b 


Of particular interest here are the theoretical predictions that dense star clust ers co uld 


h ost IMBHs formed either slowly through the successive mergers of stellar-mass BHs (Mille r fe Ham ilton 
2002h or rapidly through the runaway collisio ns of massive stars in the cluster centers 
dGurkan et al.ll2004J: iPortegies Zwart et al.ll2004j) . In addition, if IMBHs exist in star clusters, 
the inspiral of stellar-mass BHs into IMBHs i s expected to be a key source of gravitational 
waves for planned space missions (e.g., IKonstantini d is et alil2013lh Overall, finding IMBHs 
in star clusters could validate a theoretical formation channel, bolster gravitational-wave 
predictions, and test scaling relations between stellar systems and the central BHs they 
host. 


At distances of about 10 Mpc, future near-infrared telescopes will be able to measure BH 
masses as low as Mrh ~ 10 5 M & by resolving their spheres of influence in their host stellar 


systems (Do et ah 2014). While this is expected to lead to a robust inventory of IMBHs in 
the local Universe, these next-generation facilities are still many years off. Moreover, such 
sphere-of-influence studies have already been app lie d to star clusters in the Local Group, 
and all these studies are controversial ( Strader et al. l2012al . and references therein). 


Due to these on-going controversies and not wanting to wait for future facilities, we 
are pursuing an independent route, namely searching for radio_signatures of accretion onto 
putative IMBHs in extragalactic star clusters (jMaccarone 2004). Analogies with stellar- 
mass BHs suggest that an IMBH will spend more time in the hard X-ray state associated 
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with a low accretion rate, than in the soft X-ray state associated with a high accretion rate 
(reviewed by Fender &: Belloni 2012). If only a few radio observations are available, they 
are more likely to sample the steady radio emission associated with the hard X-ray state, 
as compared to the flaring radio emission associated with a transition from the hard X-ray 
state to the soft X-ray state. We thus take three approaches, as follows: 


1. Radio emission resembling that from HLX-1 in its hard X-ray state flCseh et al.l 120 15) 
could be directly detected from an individual cluster. 

2. Using the empirical fundamental-plane regression for the hard X-ray state, c onstr aints 


2003; 

Falcke et al. 

2004; 

Plotkin et al. 

2012) 


3. A conservative, semi-empirical model can be used to predict the mass of an IMBH 
that, if undergoi ng Bondi accretion in the hard X - ray state, is consist ent with the 
radio luminosity (IMaccarone fe Servillatl I2008L l2010l; IStrader et al.ll2012al) . 


In this paper, we use the Karl G. Jansky Very Large Array (VLA; 

Perlev ei 

al. 

2011 

) to 

search for radio emission from hundreds of candidate star clusters localized by 

Forbes et al. 


( 20141 ) in a mosaic of NGC 1023 with the Advanced Camera for Surveys (ACS) on the Hubble 
Space Telescope (HST). We focus on NGC 1023 because it is_one of only a few early-type 
galaxies with such an extensive ACS mosaic (lForbes et al. 2014) and, at a distance of 11.1 
Mpc, is amongst the nearest galaxies in the extensive VLA survey of K. Nyland et al. (2015, 
in preparation). At this distance 1" = 53.8 pc. 


After folding in prior information ( 

Larsen & Brodie 

2000; 

Larsen 

2001; 

Larsen & Brodie 

to 

o 

o 

to 

Brodie & Larsen 

2002; 

Morganti et al. 

2006; 

Chies-Santos et al. 

201 

3), 

Forbes et al. 


(2014) divide their candidate star clusters by age: stars older than 5 Gyr dominate the light 
of 92% of the clusters, whereas stars with ages of a few hundred Myr dominate the light of the 
remaining 8%. Those latter intermediate-age clusters may have formed during a recent gas- 
rich interaction between NGC 1023 and its companion NGC1023A, with a projected offset 
of 2(5 (8 kpc). Still, the vast majority of the clusters are dominated by an old, globular- 
clust er-like stellar popula tion. Using that pop ulation’s z-band (F 850LP) mass-to-light ratio 
from J ordan et al. (12007 ) and magnitudes from Forbes et al. (2014), the stellar masses of the 
clusters range from M* ~5x 10 4 M 0 to M* ~4x 10 6 M 0 . It is also important to underscore 
that only a handful of the candidate clusters have had their association with NGC 1023 or 
its companion confirmed spectroscopically. Liu (2 0111 ) used Chandra data to localize several 
dozen X-ray sources but did not attempt to match those sources with candidate star clusters. 
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The new VLA results are described in § [2] No individual cluster is detected, and § [3] 
explores the implications of these nondetections regarding HLX-1 analogs, X-ray detected 
clusters and semi-empirical model predictions. § Q] contains a summary and conclusions. 


2. Imaging 


The A configuration of the VLA was used to observe NGC 1023 on 2012 October 3 UT 
as part of an extensive survey for low-luminosity active galactic nuclei in early-type galaxies 
(proposal code 12B-281, PI K. Nyland). Full details will appear in K. Nyland et al. (2015, in 
preparation) but key aspects are given here. The observations were made assuming a coordi¬ 
nate equinox of 2000 and were phase referenced to the calibrator J0230+4032 at an assumed 
position of a(J2000) = 02 h 30 m 45!7108 and h(J2000) = 40°32'537068 with one-dimensional 
errors at lcr of 2 mas. The switching angle between NGC 1023 and J0230+4032 was 2?9. 
One observation of NGC 1023 was immediately preceeded and followed by observations of 
J0130+0432, leading to a net exposure time on NGC 1023 of 264 s. Every 1 s the correlator 
generated 1024 contiguous 2 MHz channels, resulting in a total bandwidth of 2.048 GHz per 
circular polarization centered at a frequency of 5.5 GHz. Observations of 3C48 were also 
acquired to set the amplitude scale to an estimated accuracy of about 3%. All but one of 
the 27 antennas provided data of acceptable quality. Given the observing strategies, the 
one-dimensional astrometric error at lcr is estimated to be 071. 


Release 4.1.0 of the Common Astronomy Software Applications (CASA) package (lMcMull in et aL 


l2007l l was used to calibrate and edit the visibility data in version 1.2.0 of an automated 
pipeline!)] The CASA task clean was used to form and deconvolve an image of the Stokes I 
emission from NGC 1023 that spanned 8000 x 07075 = 10', well matched to the HST mosaic 
flForbes et al. 2014 ). The CASA task clean was run with (1) a robustness parameter of 0.5 
to obtain the best compromise among sensitivity, spatial resolution and sidelobe suppres¬ 
sion; (2) an nterrns parameter of 2 to account for the large fractional bandwidth; and (3) the 
gridmode parameter set to “widefield” and the wprojplanes parameter set to 128 to correct 
for the effects of non-coplanar baselines. 


The CASA task widebandpbcor was then run to correct the VLA image for primary 
beam attenuation, characterized by a half-power diameter of 8'2 where the attenuation is 
50% and thus the correction is 200%. For comparison, at that location the response to a 
point source is estimated to be reduced by about 2% due to bandwidth smearing across 2 


1 https://science.nrao.edu/facilities/vla/data-processing/pipeline 
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MHz channels and by less than 1% due to time averaging over 1 s integrations^ and no 
corrections were made for these negligible effects. A square image spanning the half-power 
diameter of the primary beam on each axis was used for subsequent analysis. 

This VLA image at 5.5 GHz captured 339 star cluster candidates identified and cataloged 
by Forbes e t ah (2014). For each of those cluster s, the 2015_December 31 release of NRAO’s 
Astronomical Image Processing System (AIPS; iGreisenl 20031 ) was use d to form a cutout, 
with task subira, spanning 8" and centered on the optical position (Forbes et al. 2014). 
The ACS images were astrometrized by matching the extracted sources with CFHT archival 
catalogs with one-dimensional errors at lcr estimated to be 0"_15_(V. Pot a, 2014, priv. comm.). 
The cutout for cluster ID 39 in the nomenclature of Forbes et ah (2014) was dropped from 
further analysis becaus e it was contaminated by a 5.5-GHz counterpart to a background 
QSO ([Larsen & Brodie 2002). The radio and optical positions of this background QSO 
differed by 0"2, consistent with the error estimates quoted above for the radio and optical 
systems. ID 58 was located too close to the edge of the VLA image to form a cutout and was 
also dropped from further analysis. After this culling, 337 cutouts remain for an analysis 
of individual clusters (Figure 1). The la rms noise level amongst the cutouts has a range 
of 16-68 n Jy beam -1 , reflecting the primary beam attenuation. To minimize the risk of a 
false-positive detection of one or more individual clusters when examining an e nsemble of 
337, we adopt an individual one-tail detection threshold of 3.5cr (Wall & Jenkins 2003). 


3. Implications 


None of the 337 candidate clusters in NGC 1023 is detected above its local 3.5a level at 
5.5 GHz (Figure 1). The corresponding 3.5a radio luminosities are vL v = Lr < 4.5 —19xl0 34 
erg s -1 (Figure 2). Such a luminosity definition implicitly assumes a flat radio continuum 
spectrum up to 5.5 GHz. 


For the 311 o ld clusters, all accreting stellar-mass compact objects, whether in the hard 
X-ray state ( Strader et al. 2012bh or associ ated with a transition from the hard X-ray state 
to the soft X-ray state (ICorbel et al.l 2012|), would have radio luminosities too faint to be 
detected. For the 26 intermediate-age clusters, no other data are available for comparison. 
However, Ke p lev et ah (120141) showed that the massive star clusters in II Zw 40 at a distance 
of 10.5 Mpc have shed their enshrouding nebular emission by an age of 10 Myr. Specifically, 
the 10-Myr old cluster SSC-South in IIZw40 is not detected at 4.9 GHz and has a radio- 
luminosity upper limit of Lr < 2.1 x 10 34 erg s -1 . This strongly suggests that NGC1023’s 


2 https://science.nrao.edu/facilities/vla/docs/manuals/oss/performance/fov 
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ten-times-older clusters should lack contamination from free-free nebular emission. Thus, 
overall, we expect no contaminants for the star clusters in Figure 1, allowing us to explore 
the implications of the nondetections within the contexts of analogs of HLX-1 (§ 13. H . X-ray 
detected clusters (§ 13.2[) . and predictions of a semi-empirical model (§ 13.3p . 


3.1. HLX-1 Analogs 


14—5 


with a stellar mass of M+ ~ 10 5 6 ( 

Farrell et al. 

2009 

: Soria et al. 

2010; 

Wiersema et al. 

2010; 

Farrell et al. 

2012; 

Soria et al. 

2012; 

Farrell et al. 

2 

014 

). The stellar mass of the host 


cluster resembles those of the clusters in NGC1023, but the age of the host cluster’s light- 
dominating stars, about 20 Myr, is significantly younger. 

When transitioning from its hard to so ft X-ra y stat es, HLX-1 can achieve a luminosity 
at 7 GHz of L r ~ 3.4 x 10 36 erg s' 1 ( Webb et ah 2012). Such a level is incompatible with 
the 5.5-GHz upper limits for individual star clusters in NGC1023 (Figure 2). 

In its hard X-ray state, a combination of four measurements of HLX-1 yields a luminosity 
at 6.8 GHz of Lr ~ 1.6 x 1 0 36 er g s -1 , with no significant evidence for variability over 


several days or months ( Cseh et al.l 2015). This steady luminosity is also i ncomp atible with 
the upper limits for individual star clusters in NGC1023. As [Cseh et al.l (1201511 argue, this 
steady emission from HLX-1 is likely Doppler boosted by a factor of about five to ten, 
implying a side-on luminosity of about Lr ~ 1.6 — 3.2 x 10 35 erg s _1 . If such a side-on 
luminosity was present in any of the 337 individual clusters in NGC 1023, Figure 2 shows 
that it could escape detection in only a single cluster. 


It is important to highlight the advantages of radio searches for steady emission from 
HLX-1 analogs: its hyperluminous X-ray luminosity of Lx ~ 10 42 erg s^ 1 is _achieved only 


during a 2-3 week soft state that occurs approximately yearly (Godet et ah 20141) and can 
thus be missed if only a few X-ray observations are made. In contrast, steady radio emission 
should be present during most of the year between soft X-ray states, making it easier to 
detect even in a single radio observation, the situation for NGC 1023. Also, at the distance 
of NGC 1023, the hyperluminous X-ray luminosity of HLX-1 is below the sensitivity threshold 
of current and past X-ray all-sky monitors, blocking those discovery routes. 

Even more stringent constr aints on HLX-1 analogs in NGC 1023 come from the weighted- 
mean image stacks (e.g., Lindroos et al. 2015 1 in Figure 3 for the 311 old clusters and in 


Figure 4 for the 26 intermediate-age clusters. The stacks, which measure the mean contribu¬ 
tion to the total radio emission from each cluster population, were made using the AIPS task 
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comb. The stacks’ 3a upper limits correspond to radio luminosities of L B (old) < 2.9 x 10 33 
erg s -1 and Lr( intermediate age) < 1.2 x 10 34 erg s _1 . These stacks rule out emission from 
a population of HLX-1 analogs amongst either cluster population. 


The light-dominating stars in the clusters in NGC 1023 are at l east a factor of ten 
older than the 20-Myr-old stars in the cluster hosting HLX-1 OFarrell et a l. 2014). So one 
could speculate that if IMBHs exist in NGC 1023 clusters, those older clusters have lost 
their IMBHs due to graviational wave recoil. However, for an IMBH as massi ve as HLX-1, 
Mbr 10 4 5 Mq, ejection from the cluster is very improbable (Holley- Bockelma/nn et ajJ 


2008 b Instead, the absence of HLX-1 analogs in NGC 1023 suggests that HLX-1 is accreting 


gas somehow related to the young stars in its host cluster. An enhanced gas supply could 
arise from material remaining after the burst of star form ation 20 Myr ago, o r from material 
donated by a young stellar companion, as advocated by Miller et al, f 2014 j f 3 ! Either way, 
such an enhanced fuel supply would not be available to the IMBHs, if any exist, in the star 
clusters in NGC 1023. Alternatively, it might be that the HLX-1 phenomenon is so rare that 
a radio analog is simply not expected in NGC 1023. 


3.2. X-ray Detected Clusters 


Cluster IDs 153, 176 and 305 in t 


re nomen clature of iForbes et ajJ (120141 ) have X-ray 


detections, all of hard-state emission (Liu 20 1 lh. Using the empirical fundamental-plane 


regression for_the contracted sample of P lotk in et abj (2012), the clusters’ X-ray luminosities 
from Liu) (20111) and radio-luminosity upper limits from Figure 2 imply individual 3<r IMBH 
masses of M B h < 10 6 M 0 . These clusters have old stellar populations and estimated stellar 
masses of iff* ~ 1.3 x 10 6 M 0 or less. But a star cl uster is unlikel y to host an IMBH whose 
mass is comparable to, or exceeds, its stellar mass (Miller & Hamilton 120021 ). so these IMBH 
mass constraints are not physically interesting. Future imaging that reaches clusters with 
lower X-ray luminosities would be especially beneficial. 


3.3. Semi-Empirical Model 


Following iMaccarone &; Servillatl (120081120101) and lStrader et al.l (I2012al) . we use a semi- 
empirical model to predict the mass of a putative IMBH that, if undergoing hard-X-ray-state 
accretion in a gas-poor cluster, is consistent with the upper limit on the radio luminosity. 


3 20 Myr is approximately the main sequence lifetime of a 12 Mq star (e.g., Ford fe Bregman 20131) . 














































Specifically, we conservatively assume that the IMBH undergoes accretion at a fraction 
fb = 0.03 of the Bondi rate from a tenuous, globular-cluster-like medium with a gas density 
p = 0.2 cm -3 . Th e rati onale f or t hese assumptions and their uncertainties are discussed in 
section 3.4 of Strader et al. (2012a). These parameters yield a prediction for the hard-state 
X-ray luminosity Lx- The empirical fundamental-plane regression for the contracted sample 
of lPlotkin et al.l (120121) is then used to predict th e associated radio luminosity L\> y . (Using an 
earlier regression, Maccarone fe Serv il lat ( 20101) demonstrate that M BH oc L^ 38 (f b p)~ 0A6 .) 
A detection or upper limit to a radio luminosity thus maps to a detection or upper limit to 
an IMBH mass. 


The results of applying this semi-emipirical model to NGC 1023 are shown in Figure 
5. To avoid clutter only the best 3.5cr constraint from Figure 2 for individual clusters in 
NGC 1023 is shown. Figure 5 also conveys the 3cr constraints from the stacks of the old (Fig¬ 
ure 3) and the intermediate-age (Figure 4) clusters in NGC 1023. For comparison purposes, 


( 2008 


(20H 


in the Milkv Way Knapp et al. 

(1 

996, Pal 2, 

1GC1851, NGC 6440, 

De Riicke et al. 

(2006, NGC 6397) 

: Bash et al. 

(2008 

, M80); 

Maccarone & Servillat 


7 - - -- - -< 7 - - , f 1 I ■■■■»» ~ | I ' | ■ - ~ i 

, NGC2808); Cseh et al. ( 2Qlol. NGC 63881; Wrobel et al. ( 2011, M54) ; Lu fc Kong 


, 47Tuc, c o Gen); 


Strader et al 


(I2012al. M15, M19, M22); Ichomiuk et al.l (hohi M621: 


using distances from lHarrisI (119961) ] , as well as to the most luminous cluster in M31 (IMiller-Jones et al. 


20121 . Gl) using the distance from IVilardell et al.l (j2010h and the most lum inous cluster in 


— -X--— 

M81 

Swartz et al. 

2003 

(1994 

)• 


Mayya et al. 


20131. GC1) using the distance from Freedman et al 


From Figure 5, the best radio luminosity constraint for individual cluster candidates in 
NGC 1023 implies an IMBH mass of M B h < 3.9 x 10 5 M 0 . The constraints from the stacks of 
the old and the intermediate-age clusters in NGC 1023 correspond to mean IMBH masses of 
M BH (old) < 1.4xlO 5 M 0 and M BH (intermediate age) < 2.4 xlO 5 M 0 , respectively. However, 
recalling that the clusters’ stellar masses are M* 0.5 — 40 x 10 5 M 0 , IMBH constraints 
are only physically interesting amongst the most massive clusters with, say, about twice as 
much mass in stars as in a putative IMBH. 


Applying a z-band cut of 20.5 selects the most massive clusters, all old and with stellar 
masses A/* > 7.5 x 10 5 M 0 ( Forbes et al.1 2014). The AIPS task comb was used to form Figure 
6, the weighted-mean stack of those massive clusters. That stack’s 3cr radio luminosity is 
L B (massive) < 1.0 x 10 34 erg s _1 . Figure 5 then yields an inferred mean IMBH mass of 
M B n(massive) < 2.3 x 10 5 M 0 for the population of massive clusters. This is a physically 
interesting constraint as it suggests black-hole mass fractions of M BB (massive)/M* < 0.05 — 
0.29 for the massive star clusters. Cast somewhat differently, the clusters’ mean stellar mass 
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is 1.4 x 10 6 M 0 , implying a mean black-hole mass fraction of less than 0.16. Still, it should 
be kept in mind that the interpretation of Mbh (massive) is affected by unknowns like the 
distribution function of IMBH masses and the fraction of clusters occupied by an IMBH. 


Three massive c lust ers have been spectroscopically confirmed as being associated with 
NGC 1023 (jForbes et ah 2014) and merit special note. ID 33 and ID 269 are each as luminous 
as Gl. Assuming a uniform IMBH mass distribution and a high IMBH occupation fraction, 
ID 33 and 269 each have Mbh (massive) /M* < 0.05; for comp arison, the individual values 
inferred for Gl are ~ 0.02 from a sphere-of- influence approach (IGebhardt. Rich fe Holl20051 ). 
< 0.01 from a fundamental-plane approach ( Miller-Jones ct ah 20121 ) and < 0.02 from Figure 
5. Under the same assumptions, ID 102, an ultra-compact dwarf, has Mbh (massive)/M* < 
0.23; for comparison, Seth et ah ( 20141 ) reported an individual value of 0.15 for the ultra¬ 
compact dwarf M60-UCD1 from a sphere-of-influence approach. These comparisons suggest 
that pushing M B H(massive) lower by factors of a few could reveal analogs of Gl and M60- 
UCD1 in NGC 1023. 


Deeper VLA imaging of these massive clusters in NGC 1023 is feasible but, already, 
this stacking analysis is sensitive to the mass regime targeted for dynamical studies of_ex- 
tragalactic star clusters with next-generation near-infrared telescopes (e.g., Do et al. 2014). 
Moreover, these and additional massive clusters could be surveyed to a uniform depth in a 
future VLA mosaic that covers the full diameter, 12(6 (41 kpc ), of NGC1023’s system of 
about 500 star clusters ((Young et al. 2012; Kartha et ah 2014). In addition, the method¬ 
ology demonstrated in this work on NGC 1023 could be applie d to other nearby early-type 
galaxies with rich systems of star clu sters fe.g.. iBrodie et al.lI2014T) . Fina lly, although the 
conservative, semi-empirical model of iMaccarone fe Servillatl (120081.1201011 is appealing be¬ 
cause of its simplicity, as the photometric data impr ove o ne could consider applying more 
complex models with additional free parameters (e.g., (Sun et ah 2013). 


4. Summary and Conclusions 

We used VLA data at 5.5 GHz to search for the radiative signatures of IMBH accretion 
from 337 candidate star clusters in NGC 1023, an early-type galaxy at a distance of 11.1 Mpc. 
None of the individual clusters were detected. Similarly, only upper limits were obtained from 
weighted-mean image stacks of the 311 clusters with ages older than 5 Gyr, the 26 clusters 
with ages of a few hundred Myr, and the 20 clusters with stellar masses M* > 7.5 x 10 5 M 0 . 
We explored the implications of these data for IMBHs, if any exist in these clusters. Our 
principal findings are as follows: 
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1. The 337 candidate clusters in NGC 1023 lack radio analogs of HLX-1, a strong IMBH 
candidate in a cluster in the early-type galaxy ESO 243-49. This suggests that HLX-1 
is accreting gas related to the formation and/or presence of the 20-Myr-old stars in 
its host cluster. Alternatively, the HLX-1 phenomenon could be so rare that no radio 
analog is expected in NGC 1023. 

2. Three candidate clusters in NGC 1023 exhibit hard-state X-ray emission. From the 
empirical fundamental-plane relation, the clusters’ X-ray and radio luminosities suggest 
individual IMBH masses of Mbh < 10 6 M & . Given the clusters’ stellar masses this is 
not a physically interesting constraint. To make progress, a deeper X-ray survey that 
detects other clusters is needed. 

3. Using a formalism previously applied to Milky Way star clusters, the radio stack of the 
massive clusters in NGC 1023 implies an IMBH mass of (massive) < 2.3 x 10 5 M & . 
This physically interesting constraint can be improved with deeper VLA imaging. Also, 
more massive clusters could be added via a VLA mosaic that fully covers NGC 1023’s 
star cluster system. 


We thank the referee for a helpful report^Dr. Vince Pota for estimating the astrometric 
accuracy of the HST positions reported in Forbes et ah (12014]) and Dr. Craig Walker for 
advice on image stacking. JCAMJ is the recipient of an Australian Research Council Future 
Fellowship (FT140101082). JMW is grateful to the NRAO Director, Dr. Tony Beasley, for 
approving the leave of absence that facilitated this research. 
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CONT: NGC 102 5499.000 MHz NGC1023-IND.SUBIM.1 



Fig. 1.— VLA cutout of the Stokes / emission at 5.5 GHz centered on the optical position 
of a candidate star cluster in NGC 1023. The hatched ellipse in the north-east corner shows 
the beam dimensions at FWHM of Cf/46 (25 pc) x Cf/35 (19 pc) with an elongation position 
angle of 61°. Contours are at -6, -4, -2, 2, 4, and 6 times the la rms noise shown in the 
legend in units of Jy beam -1 . Linearly-spaced contours are optimal for conveying noise levels 
at a glance. Negative contours are dashed and positive ones are solid. The central circle 
of diameter O'.'6 (32 pc) shows the cluster’s optical positional uncertainty at 95% confidence 
and is labelled with its ID from Forbes et al. (2014). The VLA photometry seeks evidence 
for the accretion signature of a point-like IMBH in the cluster’s center. The cluster is not 
detected above the 3.5a level. Figures 1.1 — 1.337 are available in the online version of the 
Journal. Some cutouts show additional clusters offset from the central cluster. The 337 
clusters have effective optical diameters of 0.8 — 40 pc (Forbes et al. 2014) and their images 
are ordered by increasing cluster diameter. 
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Individual Candidate Cluster 

Fig. 2.— Upper limits to the radio luminosities at 5.5 GHz, L R , of 337 individual candidate 
clusters in NGC 1023. The abscissa conveys the ordering of the clusters, 1 through 337, and 
corresponds to the ordering of Figures 1.1 - 1.337. 
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Fig. 3.— Weighted-mean stack of the VLA images of the 311 old clusters in Figure 1. The 
stacked image has an rms noise of 1.2 /i.Jy beam -1 (lcr). The hatched ellipse, contouring 
scheme and circle are the same as for Figure 1. The stack measures the mean contribution 
to the total radio emission from the old clusters. No emission is detected above 3 a = 
3.6 /iJy beam -1 . 
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Fig. 4.— Weighted-mean stack of the VLA images of the 26 intermediate-age clusters in 
Figure 1. The stacked image has an rms noise of 4.9 /xJy beam -1 (la). The hatched ellipse, 
contouring scheme and circle are the same as for Figure 1. The stack measures the mean 
contribution to the total radio emission from the intermediate-age clusters. No emission is 
detected above 3a = 14.7 /i.Jy beam -1 . 
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Fig. 5.— IMBH masses predicted from a semi-empirical model for star clusters at various 
distances. If a putative IMBH is undergoing hard-X-ray-state accretion in a gas-poor cluster, 
then a detection of, or upper limit to, its radio luminosity yields a value for, or upper limit 
to, its IMBH mass. To visually convey the model’s mass sensitivity with distance, the dotted 
line shows the trend for a fiducial frequency of 6 GHz and a fiducial flux density of 30 /iJy. 
Upper limits achieved for NGC 1023 are shown in color. For context, upper limits are shown, 
in black, for clusters in the Milky Way, M31 and M81 (see the text). 
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Fig. 6.— Weighted-mean stack of the VLA images of the 20 most massive clusters in Figure 1. 
The stacked image has an rms noise of 4.3 /xJy beam -1 (lcr). The hatched ellipse, contouring 
scheme and circle are the same as for Figure 1. The stack measures the mean contribution 
to the total radio emission from the most massive clusters. No emission is detected above 
3cr = 12.9 /i.Jy beam -1 . 




